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ABSTRACT 

We present new BVI ground-based photometry and VI space-based photometry for 
the globular cluster M92 (NGC 6341) and examine luminosity functions in B, V, and 
I containing over 50,000 stars ranging from the tip of the red giant branch to several 
magnitudes below the main sequence turn off. Once corrected for completeness, the 
observed luminosity functions agree very well with theoretical models and do not show 
stellar excesses in any region of the luminosity function. Using reduced-^ 2 fitting, the new 
M92 luminosity function is shown to be an excellent match to the previously published 
luminosity function for M30. These points combine to establish that the "subgiant 
excess" found in previously published luminosity functions of Galactic globular clusters 
are due to deficiencies in the stellar models used at that time. Using up to date stellar 
models results in good agreement between observations and theory. 

Several statistical methods are presented to best determine the age of M92. These 
methods prove to be insensitive to the exact choice of metallicity within the published 
range. Using [Fe/H]=— 2.17 to match recent studies we find an age of 14.2 ± 1.2 Gyr for 
the cluster. 



Subject headings: globular clusters: individual(M92), stars: distances, stars: evolution 



1. Introduction 

Globular clusters are ideal locations to test 
stellar evolutionary models due to their single- 
age single-metallicity nature. Previously, 
work has focussed in great depth on analy- 
sis of the color-magnit ude diagram (CMP) of 
cluste rs, for example in Renzini k. Fusi Pecci I 
(|l988l ). While studies of the CMD can re- 
veal a great deal about stellar evolution, the 
luminosity function (LF) of the cluster is espe- 
cially powerful for determining the timescales 
of stellar evolution. Below the main-sequence 
turnoff (MSTO), the LF is primarily a reflec- 
tion of the initial mass function modulated by 
dynamical mass segregation effects at the low- 
est masses. However, above the MSTO, the 
LF reveals the progress of the hydrogen burn- 
ing shell through the star and can even give 
hints as to the internal structure of the star. 
Indeed, the enhancement in the LF known as 
the LF 'bump' marks the hydrogen-burning 
shell's transition from the region previously 
mixed by convecti on into unmixed stellar ma- 



terial (|lben 1119681 ) 



Only with the advent of large-format CCD 
cameras has it become possible to obtain pre- 
cise photometry of a large numbers of stars 
in all phases of stellar evolution. Previous 
studies have been largely limited to either 
low precision photographic measurements or 
small spatial coverage and limited samples of 
stars. With deta iled LF s , une xpected results 



have appeared. Bolte ( 19941 ) found an ex- 



cess of stars on the subgiant branch (SGB) 
for th e low metallicity, [Fe/H]=— 2.12 ([Harris 
1996), cluster M30 (NGC 7099). The sub 



gia nt excess in M30 was further confirmed 
by ISandquist et al. I ( 1999h using a higher- 
quality data set. This SGB excess is the 
expected result if weakly interacting massive 
particle (WI MP) energy transport is im por- 
tant in stars ([Faulkner fc Swenson 1993). Rc- 
cent LFs of the more m etal rich clusters M5 
(ISandquist et al. Ill996l ) and M3 (|Rood et al. 
1999) with [Fe/H]=-1.27 and -1.57 respec 



poor clusters. However, the models used in 
the comparison do not include diffusion, which 
is a standard part of modern stellar evolution 
codes and the statistical basis for the SGB ex- 
cess is not explicitly defined. 

The luminosity function of M92 (NGC 
6341) is an ideal test for stellar modeling 
codes . Its metallicity, [Fe/H]=— 2.27 ( Harris I 



19961 ), places it in the same abundance range 



as M30, which was earlier found to have a SGB 
excess. It is a large cluster, making it easy to 
measure significant numbers of stars. Its lo- 
cation, f ar from the g alactic plane (b = 34.86 
degrees (Harris II1996IU . minimizes the signif- 
icance of field-star contamination. Finally, it 
is a fairly well studied cluster with accurately 
determined distance modulus and metallicity 
which simplifies comparisons of the observed 
LF to models. Fina lly, previous LF stud ies of 



the cluster such as iPiotto et al. I (|1997l ) and 



Lee et al. 



(|2003l ) have not looked at the LF 
along the RGB in order to examine the pre- 
viously found excess present in M30, instead 
they have concen trated on the lower MS. The 
Lee et al. ( 20031 ) study further supports us- 



ing M92 as a test for stellar modeling codes 
since examination of the mass function in that 
work suggests that the cluster has not been 
strongly affected by tidal shocks resulting in 
a pure sample of cluster stars. 

In the context of the above points, we 
herein present our analysis of a comprehen- 
sive luminosity function study of M92 us- 
ing ground-based and space-based observa- 
tions. The next section describes these ob- 
servations, while Section 3 discusses the re- 
duction of the data. Sections 4 and 5 present 
the color-magnitude diagrams and luminosity 
functions. Section 6 details the stellar evo- 
lution model and discusses details in fitting 
the theoretical models to the observed LFs. 
Section 7 compares the new M92 LF to pre- 
vious LFs of M92 and M30 and Section 8 
presents conclusions about the fits and the 
general state of LF modeling. 



tively do not show the SGB excess, suggesting 
that it may only be a characteristic of metal- 



1 



2. Observations 

The ground-based images were collected 
from 21-27 April 2003 using the Hiltner 
2.4m telescope at MDM Observatory on Kitt 
Peak. All frames were taken using the Echellc 
2048x2048 camera with a pixel scale of 0.28 
arsec/pixel and a total field of view of 9.56 
arcminutes. To cover a significant fraction 
of the cluster area, short and long exposures 
were tiled over a 27x27 arcminute area. For 
the short images, exposure times were cho- 
sen to ensure that stars at the tip of the red 
giant branch (RGB) were not saturated, typi- 
cally resulting in 5-15 second exposures. Two 
images in each filter were taken at each grid 
position for the short images. The exposure 
lengths for the longer exposures were in the 
range of 60-120 seconds depending on the fil- 
ter and 3 exposures were taken in each filter 
at each grid position. All of the images were 
observed under photometric conditions with 
average seeing of approximately 0.8 arcsec- 
onds. 

The Hubble Space Telescope images were 
taken from the STSCI archives and were col- 
lected on 27 August, 2002, using the Ad- 
vanced Camera for Surveys. The Hubble 
Space Telescope images were obtained with 
the Advanced Camera for Surveys on 27 Au- 
gust 2002 as part of HST program GO-9453. 
The exposures cover the central 3-arcminutes 
of the cluster in the F606W and F814W filters, 
with exposure lengths of 0.5, 5, and 90 sec- 
onds and 0.5, 6, and 100 seconds respectively 
in the two filters. Because these observations 
were taken shortly after ACS was installed, 
the charge transfer efficiency corrections are 
neglig ible (|Brown et al. Il2005h . 



3. Data Reduction 
3.1. Image Processing 

The ground-based images of M92 were pro- 
cessed using standard IRAF techniques uti- 
lizing evening and dawn twilight flats. The 
archival Hubble images were processed us- 
ing the STSCI on-the-fly-reprocessing system 



from the Multi-mission Archive at Space Tele- 
scope (MAST). The drizzled images, multi- 
plied by the image exposure time, were uti- 
lized for the photometry. 

3.2. Photometry 

The ground-based sample 

Profile fitting photometry was performed 
on the 179 image frames using the DAOPHOT 
and ALLSTA R programs develop ed by Pe- 

In each 
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tcr Stetson ([Stetson 111987 . 
frame, approximately 100 bright uncrowded 
stars were chosen to determine the point- 
spread function (PSF) and its variation about 
the frame. While the PSF stars were specif- 
ically chosen to be outside the most crowded 
regions of each frame, they generally did sam- 
ple a large fraction of the frame and therefore 
map the PSF variations well. Several scripts 
were used to automate two passes through 
DAOPHOT 's FIND routine and ALLSTAR to 
generate a photometry list for each image. 

Photometry from the individual frames was 
filtered to remove any detection with a mea- 
sured error greater than 0.1 mag. Aperture 
corrections were calculated using the bright- 
est uncrowded stars on each frame. These 
were used to search for a spatial dependence of 
the aperture correction, but none was found. 
DAOMASTER and DAOMATCH were then 
used to combine the individual photometry 
files into one master file for each filter, re- 
quiring that a star be detected in at least two 
frames in each filter. These master files were 
filtered to contain only stars whose frame-to- 
frame magnitude variation was less than 0.1 
magnitude. The B, V, and I master files were 
then matched using DAOMASTER and DAO- 
MATCH requiring a star to be detected in all 
three filters to be included in the final catalog. 

The HST sample 

The space-based images were processed us- 
ing the same methods as the ground based 
data, although the process was simplified by 
the fact that there were only 6 images. While 
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the same criteria for rejecting stars based on 
photometric errors was used, the matching cri- 
terion was relaxed to require that stars only 
be detected once in each filter. 

3.3. Calibration 

The instrumental magnitudes were brought 
onto the standard system using P. Stetson's 
photo metric standards for M92 (jStetson 
2000). Using stars brighter than 17th mag- 
nitude in V, we were able to find over 600 
stars with B, V, and I magnitudes in com- 
mon between the Stetson standards and our 
ground-based data. The best photometric so- 
lution to bring the data to the Stetson system 
was found to depend on color only to first- 
order. The transformation equations were 
determined to be 

B = b + 2.685 + 0.0537 (B - V) 

V = v + 1.530 - 0.0413 (B - V) 

/ = i + 0.223 + 0.0120 {V -I) 

The residuals from the fit between the ob- 
served stars and the Stetson standards are 
shown in Figure [TJ In all cases, the distribu- 
tion of the residuals is exactly the distribution 
expected from the photometric errors. 

The HST data posed an interesting prob- 
lem in calibration because of the small field 
of view in the crowded core of the cluster. 
Due to confusion between stars in the core 
of the cluster, there were no Stetson stan- 
dard stars in the field of view of the HST im- 
ages. It was possible, however, to calibrate 
the HST photometry using our ground-based 
observations as secondary standards. Using 
the same method as in the calibration of the 
ground-based data, approximately 400 stars 
in common were compared. These stars all lie 
in the outer 0.5 arcminutes of the HST im- 
ages. In order to simplify the interpretation 
of the luminosity functions, which are based 
on BVI magnitudes, the instrumental magni- 
tudes from the Hubble data were transformed 
directly into V and I using the following rela- 
tionships 

V = 606 mst +0.4881+0.0695 (V-I)+0.1510 {V-lf 



I = 814 inst -0.3850-0.0521 (U-/)+0.0306 {V-lf 

Figure [2] shows the residuals between the HST 
photometry and the ground based photome- 
try. The distribution matches expectations 
based on the photometric errors in the two 
data sets. 

4. The Color-Magnitude Diagrams 

In total, 34,242 stars in the ground-based 
sample and 41,205 stars in the HST sample 
were measured. The HST sample is larger due 
to the decreased confusion in the core which 
allowed a much larger number of MS stars to 
be detected. V versus (B — V) and (V — I) 
color-magnitude diagrams from the terrestrial 
sample are shown in Figure [3] The V ver- 
sus (V — I) color-magnitude diagram from the 
HST sample is shown in Figured] 

5. The Luminosity Functions 
5.1. Completeness 

The critical task in generating a luminos- 
ity function is determining the completeness of 
the photometric sample. Generally, the com- 
pleteness is a function of both position relative 
to the cluster, due to the confusion in crowded 
regions, and magnitude. An accurate determi- 
nation of the completeness is found through 
extensive use of artificial star tests. These 
tests were performed by adding artificial stars 
to each image and then remeasuring photom- 
etry for the complete field to determine what 
fraction of the added artificial stars are recov- 
ered. Errors in the eventual LFs are governed 
by Poisson statistics, so large numbers of in- 
put and recovered artificial stars are needed. 
However, adding a large number of stars to 
any particular image will change the crowd- 
ing and therefore the completeness. In order 
to balance this effect, many runs through the 
artificial star routine are used with relatively 
small numbers of stars added in each run. 

The general method used for the artificial 
star tests was as follows: a master list of artifi- 
cial stars with random x and y pixel positions 
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Fig. 1. — Residual distributions for the stan- 
dard stars. The second column of plots shows 
the error distribution displayed with fitted 
gaussians with sigma 0.023, 0.028, and 0.033 
for B, V, and I. These sigmas are as expected 
for simple error propagation given the photo- 
metric errors. 




Fig. 2. — Residual distributions for the stan- 
dardized HST sample. The residual distribu- 
tion is shown overplotted with gaussians with 
standard deviations of 0.05 in V and 0.06 in I. 




Fig. 3.— The V - (B - V) and V - (V - 
I) color-magnitude diagrams for the ground- 
based sample. 
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Fig. 4. — The V - (V — I) color-magnitude di- 
agram for the HST sample. This CMD covers 
an area of approximately 3x3 arcminutes and 
contains over 41,000 stars. 
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covering the entire area imaged was created. 
The magnitude for the artificial stars was de- 
termined by generating random numbers cov- 
ering the range of instrumental B magnitudes 
found during the photometry. The assigned 
B magnitudes were then used with the CMD 
ridge lines to assign V and I magnitudes to 
each artificial star. The result is a list of ar- 
tificial stars distributed randomly across the 
imaged area with magnitudes and colors that 
match the observed cluster stars. Each star 
from the master list was placed at the appro- 
priate position with the appropriate instru- 
mental magnitude in each frame that it could 
have appeared. The image frames with added 
artificial stars were then put back through 
the same DAOPHOT and ALLSTAR pipeline 
used for the initial photometry including the 
same error clipping and required number of 
detections during matching. Additionally, the 
input and recovered magnitudes for the artifi- 
cial stars were required to be within 0.1 mag 
in order to remove the possibility of real stars 
or blends being confused with the artificial 
stars. Figure [5] shows the difference between 
the input and recovered artificial star magni- 
tudes as a function of V magnitude. To in- 
vestigate any magnitude bias in the recovered 
artificial stars, the stars were divided into bins 
two magnitudes wide and the average differ- 
ence between the input and recovered mag- 
nitude was determined. For all bins brighter 
than 22th magnitude, the bias had an abso- 
lute value of 0.002 magnitudes or less. The 
final bin had bias of 0.012. However, the last 
bin is beyond the completeness limit and has 
an order of magnitude less stars than any of 
the other bins. Based on this analysis, there 
is no significant bias in the artificial star mag- 
nitudes. In all, 24 sets of artificial stars, tak- 
ing approximately 36 cpu- hours per set, were 
added to the images for a total of over 91,000 
input stars and 59,000 recovered stars. 

Due to the relatively small number of im- 
ages in the space-based data, it was possible 
to use a larger number of artificial star runs in 
a reasonable amount of time. Using the same 
procedure, over 267,000 artificial stars were 



added and over 203,000 were recovered from 
the images after processing 90 sets of artificial 
stars. The completeness of the datasets as a 
function of magnitude is shown in Figure [S] 
The lower completeness at the bright end is 
due to bright stars saturating in some of the 
images. 

5.2. Computing the Luminosity Func- 
tion 

There are two main steps in transforming 
the photometric data into the luminosity func- 
tion. First, field stars, horizontal branch (HB) 
stars, and asymptotic giant branch (AGB) 
stars must be removed to produce a pure sam- 
ple of stars on the MS and RGB. The field 
stars, HB stars, and AGB stars were removed 
by clipping all stars more than 3-cr in color 
from the fiducial reference line. To find the 
fiducial line, a box was run across the CMD 
in color at different magnitudes. The set of 
points in color that maximize the population 
in the box as a function of magnitude define 
the fiducial line. This method will not remove 
any field stars that lay along the MS and RGB, 
but the density of field stars is low enough that 
the resulting LF is not significantly affected. 
This is confirmed by the B esancon stellar pop- 



ulation synthesis models (jRobin et al. 1 12003) 



which predict approximately ten field stars 
along the upper half of the RGB. The pho- 
tometric clipping method is ambiguous at the 
top of the RGB where the AGB and RGB 
nearly merge. In this region, extraneous AGB 
stars were removed manually. 

Completeness corrections were accom- 
plished by weighting each star individually. 
A grid of completeness as a function of ra- 
dius from the cluster center and magnitude 
was constructed and then used to assign a 
weight of completeness(m, r) _1 to each star. 
The LF was then constructed by summing the 
weights of stars in each magnitude bin. The 
completeness limit was defined as the magni- 
tude where the average weight per star in the 
bin was equal to 2, roughly the point where 
the data is 50% complete. The completeness 



5 



0.05 

1 

1 

*. o 



T. r 




Fig. 5. — The difference between input and 
recovered artificial star magnitudes for the 
59,000 recovered artificial stars as a function 
of V. The average bias is insignificant at - 
0.0005 magnitudes in V and similar amounts 
in B and I. 




Fig. 6. — The completeness of the ground- 
based (solid line) and space-based (dashed 
line) data as a function of V magnitude. The 
completeness limits quoted for the LFs are 
slightly different than may be expected from 
these lines since they are based on complete- 
ness as sampled by the positions of the cluster 
stars rather the random positions of the arti- 
ficial stars. 



in the V band as a function of magnitude is 
shown in Figure 6. The error bars on the 
luminosity function are Poisson errors given 
by 



a(Ni) = Niy/l/m + l/ai + l/bi 

where Ni is the summed weights of stars in 
the magnitude bin, is the number of stars 
in the magnitude bin, ai is the number of re- 
covered artificial stars in the magnitude bin, 
and bi is the number of input artificial stars 
in the magnitude bin. The number of stars in 
each bin is the dominant source of uncertainty 
rather than the number of input or recovered 
artificial stars. 

This weighting method revealed one short- 
coming of the ground-based data, namely the 
fact that, due to the high stellar density, no 
stars fainter than the SGB were detected in 
the central two arcminutes of the cluster. This 
introduces a skew into the LF since the area 
being sampled at each magnitude is differ- 
ent. This could be corrected by defining a 
cluster profile based on brighter, more com- 
plete, stars and then correcting the numbers 
of fainter, less complete, stars. This method 
has been used in previous papers in the liter- 
ature. However, any correction would assume 
that the cluster has the same profile for high 
and low mass stars and would completely dis- 
regard any mass seg regation in the clus ter, as 
was found for M3 bv lRood et al~l (|l999l ). This 
method could also cause problems since errors 
in the numbers of bright stars, with poor er- 
ror statistics, would be propagated through 
the entire LF. Also, because the size of the 
central hole depends very strongly on magni- 
tude, extreme care would be required to avoid 
artificially generating an excess or deficiency 
at points in the LF. To simplify the analy- 
sis of the LF, the central area of the ground 
based data, containing approximately 11,000 
stars, was removed. The resulting complete- 
ness limits are B=22.1, V=21.3, and 1=20.5. 
The V luminosity functions from the ground, 
HST, and combined data sets are shown in 
Figure [7] and are representative of the LFs in 
the other bands. The ground based LFs are 
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listed in Table 1 in the appendix. 

The HST sample displayed reasonable com- 
pleteness across the entire field, unlike the 
ground based sample, so the entire catalog was 
used to create the luminosity functions using 
the procedure outlined above. Due to the rel- 
atively small field of view of the ACS, very few 
stars near the tip of the RGB were detected 
resulting in a very noisy LF. The completeness 
limits arc V=21.1 and 1=20.1, limited largely 
by the short exposure times. The LFs from 
the HST sample are listed in Table 2 in the 
appendix. 

Combining the HST data from the cluster's 
central region with the ground-based pho- 
tometry of the outer region, it was possible 
to create a composite LF covering a full 27 
arcminute square area of the cluster. The 
composite LF contains data for over 54,400 
weighted stars. The combined LFs are listed 
in Table 3 in the appendix. The complete- 
ness limits for this sample are taken to be the 
completeness limits of the HST sample. 

6. Theoretical LF Models and Fitting 

To produce the theoretical LFs for compar- 
ison to the observed data, stellar models were 
created using the D artmouth Stellar Evolu- 
tion Code (DSEP) dChabover et al. I l200lh. 



These models used the lVandenBerg and Clem 
(|2003l ) color transformations. A catalog of 
theoretical luminosity functions was created 
with ages stepped by 0.1 Gyr from 11.6 Gyr 
to 16.5 Gyr. The models use scaled solar com- 
positions, while globular clusters are typically 
enhanced in their a-element (O, Mg, Si, S, and 
Ca) a bund ances. As noted by IChieffi et al. I 
(|l99ll ) and IChabover et~aLl (jl992j), scaled so 
lar composition models are nearly identical 
to a-element enhanced models, provided one 
modifies the relationship between [Fe/H] and 
the heavy element mass fraction Z . The mod- 
ification assumed [a/Fe] = +0.40. Metallici- 
ties of [Fe/H]= -2.11, -2.14, -2.17, -2.20, 
—2.23, —2.27, and —2.31 were chosen to 
sample the publishe d range of —2.24 ± .0 8 
(IZinn fc West I Il984h . -2.27 (|Harris I fl~996h . 




Fig. 7. — From top to bottom: the ground- 
based V LF, the HST V LF, and the combined 
V LF. In all three LFs, the error bars are the 
1-er Poisson errors. The LFs in B and I look 
very similar and are of equal quality. 



7 



and -2.38 ± 0.07 (jKraft fe Ivans Il2003h . The 
majority of the comparisons between the ob- 
served and theoretical LFs were completed 
using the [Fe/H]= —2.17 models since they 
are the closest metall icity match to the recent 
Kraf t fc Ivans 1 (|2003h measurements after be- 
ing modified for the scaled-solar assumption. 
All of the statistical tests proved to be insen- 
sitive to metallicity in the studied range. 

One of main parameters in fitting the the- 
oretical LFs to the observed ones is distance 
modulus. To ensure that this represents the 
absolute distance modulus, the observed LFs 
must be corrected for extinction. A stan- 
dard extinction-color excess relation was used, 
A(V) =3.1 E(B - V), with the published 
value of the c olor excess being 0.02 for M92 
([Harris 1 11996). Extinction in the I band was 
calculated assuming E(V-I) = 1 .25 E(B-V) 
([Dean. Warren, fc Cousins I ll978h . All extinc- 
tions wer e computed on t he B, V, and I mag- 



nitudes. ISirianni et al. I (|2005l ) suggest that 



reddening and extinction should be consid- 
ered in the native 606W and 814W filters for 
the space-based observations due to the slight 
differences with the standard V and I filters. 
Due to the low reddening to M92, the differ- 
ence in methodology corresponds to a differ- 
ence of 0.008 magnitudes in the V filter and 
0.001 magnitudes in the I filter. 

Three separate statistical methods were 
used to determine the best match between 
the observed and theoretical LFs. First, a 



the tip of the RGB to one magnitude be- 
low the MSTO, to V = 18.9. This covers a 
mass range from 0.73 to 0.77 Mq and there- 
fore tests the relative numbers of stars on the 
RGB and MS rather than being sensitive to 
the cluster IMF. Standard probability tables 
for the reduced-x 2 statistic assume Gaussian 
errors while the errors in the LF come from 
Poisson counting. As a result, we employ a 
Monte Carlo simulation to estimate the sta- 
tistical significance of our LF comparisons. In 
the simulation, 10 6 realizations of the M92 
LF were generated using the quoted errors. 
The reduced-x 2 values between the original 
LF and the realizations were then computed 
to generate a probability table. Reduced-x 2 
values of 0.9 or less are found 75% of the time 
showing a good match between the observed 
LF and the theoretical LFs. A typical best 
match between the theoretical and observed 
ground-based LFs is shown in Figure [5] 

6.1. The K-S Test 

In order to get a better constraint on the 
data, the observed and theoretical data were 
compared using the Kolmogorov-Smirnov (K- 
S) test. The K-S test maps the maximum 
difference between the theoretical cumulative 
LF and the observed cumulative LF and is 
therefore somewhat immune from the poten- 
tial blending of stars between bins by photo- 
metric errors in the conventional LF. 

Xq complete the K-S test, a subsection of 



reduced-x 2 fitting method based on| Bevington fc Robinfio/i lu minosity f unct ions from approximately 
(1992) was used to find the best match be- 
tween the theoretical and observed LFs in 
all three filters simulatinously. During the 



reduced-x minimization, the absolute dis- 
tance modulus was allowed to vary while the 
normalization was set by the total number of 
stars brighter than the completeness limit in 
the LF. Regardless of the theoretical model 
used, the reduced \ 2 values were similar with 
values around 0.9 with a wide range of dis- 
tance moduli (±0.3 magnitudes) and ages 
(±1.2 Gyr) producing seemingly reasonable 
fits. The reduced-x 2 fits covered a range from 



0.5 mag fainter than the tip down to the 
MSTO was defined. The cumulative LF from 
this region is then compared to the cumula- 
tive theoretical LF while stepping through the 
absolute distance modulus. Ages from 11.6 
Gyr to 16.5 Gyr were examined in 100 Myr 
increments. The K-S probability statistic was 
then used to find the best match. The error 
in age and absolute distance modulus was as- 
sumed to be given by the error ellipse traced 
by the 50% contour line. There was no signif- 
icant difference between the different metal- 
licity models so comparisons have been made 
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using the model closest to the lKraft fc Ivans" 
(|2003l ) metallicity, the [Fe/H]= -2.17 model. 
This results in good agreement between the 
observed LFs and the model giving an age of 
14.2 ± 1.2 Gyr. This result is shown in Figure 

m 

7. Comparisons 

7.1. Previous M92 Luminosity Func- 
tions 



LeeetaL 



(|2003l ) published a luminos- 
ity function for stars fainter than than the 
MSTO, which serves as a good check for the 
luminosity function presented in this paper on 
the low mass end. The shape of the Lee et al. 
luminosity function agrees very well with our 
luminosity function over the range of V=18.5 
to V=20.9 giving a reduced x 2 value of 0.52 
once the two LFs are normalized. The match 
between the LFs is shown in Figure |TD] 

7.2. The M30 Luminosity Function 



Pre vious work such as iSandauist et al 
(|l999h have shown that the LF of M30 has an 
excess of stars on the subgiant branch when 
compared with stellar models. It is not clear 
if this behavior is particular to M30 or if it 
is a general characteristic of low-metallicity 
globular clusters. To compare the LFs from 
M30 and M92, the reduced-x 2 between the 
two cluster LFs ignoring the small metallicity 
difference was calculated. To find the mini- 
mum reduced-x 2 , the distance modulus and 
normalization of M30 were allowed to vary. 
Assuming a distance modulus of 14.64 for 
M92 the best fits were found using a distance 
modulus of 14.92 for the V LF and a distance 
modulus of 14.82 for the I LF. This matches 
well with the distance modulus of 14.87 ±0.12 
found bv lSandauist et al~ (1999) assuming a 
reddening of E(V — I) = 0.06. The resulting 
reduced-x 2 values are 0.26 in V and 0.36 in I 
and suggest that the LFs for M30 and M92 are 
very similar. Given the excellent agreement 
between the theoretical LFs and the M92 LFs 
and given the extreme similarity between the 




Fig. 8.— The theoretical 13.9 Gyr [Fe/H] = 
— 2.17 LF, the solid line, plotted with the ob- 
served V LF from the ground-based sample, 
the error bars. The reduced-x 2 is 0.96. The 
B and I fits have equal quality. 




Fig. 9.— The K-S test results for [Fe/H] = 
—2.17 models. The color scale gives the con- 
fidence as a percentage from the K-S test, the 
contours are drawn at the 50% level. 
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M92 LFs and previous M30 LFs, there seems 
to be no problem with current stellar evolu- 
tion models in contrast to previous results of 
ISandauist et al. I (|l999h . The match between 
the M30 and M92 V LFs are shown in Figure 
[HI and is also representative of the quality of 
the fit in the I band. 

1 19991 ) 



Additionally, the iSandquist et al 



LF was directly compared to models with 
[Fe/H]=-2.02 and -2.42 to bracket the nom- 
inal metallicity. Allowing for a slight uncer- 
tainty in the distance modulus and normaliz- 
ing the theoretical LFs to match the observed 
LFs, reduced-x 2 values of 1.16 and 0.95 were 
found for the two models. Using the generated 
probability table based on Poisson errors, we 
find that a reduced-^ 2 value less than 0.95 is 
found 53% of the time and a value of 1.16 is 
found 25% of the time. From this, we can con- 



clude that the M30 LFs of ISandauist et al. 



(jl999l) are well fit by DSEP models. 

Adding further support to this argument 
is a direct co mparison of the luminos ity func- 
tions used in lSandquist et al~l (jl999h to cur- 
rent DSEP luminosity functions. Figure [T2l 
compares the standard DSEP model, includ- 
ing diffusion and the most recent reaction 
rate for 14 A + P -^ 15 O + 7, the slowe st 



rate in the CNO cycle, from iFormicola et al 



(|2004l ) to models with th e old and new rates 
witho ut diffusion and IVandenBerg et al. 
(|2006f ) models. All of the LFs are at the 
same metallicity and a re 14 Gyr old. The 



Sandquist et al. ( 19991) paper utilize d 



(2006) 



lier versions of the IVandenBerg et al. 
models. F rom the Figure, i t is apparent 



that the IVandenBerg et al. 
are similar to the DSEP 
els. Compared to the st a ndard DSEP model, 



( 20061 ) models 
no-diffusion mod- 



(2006) models have 



the IVandenBerg et al. 
9% fewer stars along the RGB. The differ- 
ence explains the deviation found between 
the observed LF and the theoretical LFs in 



Sandquist et al. I (|l999h 



Th is underscores a point raised in lGallart et al. 
(|2005l ) that a variety of different model li- 
braries should be utilized before assuming 




Fig. 10.— The lLee et al. I (|2003D LF and our 



LF in the V band. Our LF has 1-sigma error 
bars while the Lee LF is plotted without error 
bars for clarity. 




Fig. 11.— The ISandquist et al. I (|l999h M30 



LF and our LF with l-er error bars in the V 
band. For clarity, the error bars have been 
omitted from the M30 LF. 
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Fig. 12.— 14 Gyr DSEP LFs with and 
without current reaction rates and diffu- 
sion compared to IVandenBerg et a"l~ (.2006) 
LFs. When normalized to the MS, the 



VandenBerg et al. ( 20061 ) models have 9% 
fewer stars along the RGB. This decrease 
results in the excess of stars found in 
Sandquist et al. I (|l999h . 



some unaccounted phenomenon is responsi- 
ble for deviations between observed and the- 
oretical LFs. It also gives a great amount of 
confidence in stellar evolution models which 
use up to date physics, such as DSEP. Fur- 
thermore, the comparison between DSEP and 
VandenBerg et alH ( 20061 ) LFs and the suc- 
cesses and failures in matching observed lumi- 
nosity functions of low-metallicity GCs sup- 
ports the inclusion of gravitational settling 
and microscopic diffusion in stellar evolution 
models and the implied reduction in globular 
cluster ages of order 1 Gyr. 

7.3. The RGB Bump 

The bump is caused by an increase in the 
hydrogen content of the material fed into 
the hydrogen-burning shell in stars along the 
RGB. This increase comes when the shell 
passes throu gh the form er base of the con- 



vection zone (jlben Ill968[ ) . While it is difficult 



to make a concrete comparison due to the 
large error bars in the LF along the RGB, 
after examining the best matches it appears 
that M92's bump is approximately 0.4 mag 
fainter than the theoretical b ump. This agrees 
witfiprevious results from iFusi Pecci et al. 
(J) which found a 0.4 mag difference for 



(199C 



a sample of 13 clusters. Two possible ex- 
planations for the difference exist: either 
the hydrogen-burning shell has moved more 
rapidly than expected or convection has pen- 
etrated farther into the star than expected, 
perhaps through a mechanism such as con- 
vective overshoot. Addition ally, recent results 
( Biork and Chabover 20061 ) suggest that the 
location of the bump is extremely depen- 
dent on metallicity and composition with 0.2 
magnitude or more changes being possible 
within reasonable ranges of metallicity and 
a-enhancement . 



8. Conclusions 



1. This paper finds good agreement be- 
tween theoretical LFs and observations in 
terms of relative stars counts along the MS, 
SGB, and RGB in the M92 LF. While the 
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LFs are noisy in the RGB region, the theo- 
retical LF matches the observed LF within 
the errorbars. The result is further strength- 
ened by the use of two independant data sets: 
the ground based observations and the HST 
observations. 

2. Contrary to previously published results, 
we find no "subgiant excess" in comparisons 
of the M92 LF to DSEP LF s, the M92 LF 
to the lSandauist et al. I (| 19991) of M30, or the 



M30 LF directly to the DSEP LFs. It appears 
that the "subgiant excess" found in previous 
works is due to an erroneous under-prediction 
of the models rather than a real excess of 
stars in the cluster. Direct compariso n of 
DSEP LFs to IVandenBerg et al. I (|2006l ) LFs 



reveal that gravitational settling and micro- 
scopic diffusion all combine to better fit ob- 
servations. Thus, observations of metal-poor 
GCs support the inclusion of these effects in 
stellar evolution models and their implied re- 
duction in globular cluster ages. 

3. Utilizing two separate methods to com- 
pare the theoretical and observed LFs it is 
possible to constrain the age of the cluster to 
14.2±1.2 Gyr with an absolute distance modu- 
lus of 14. 60 ±0.09. In all cases, the comparison 
models included diffusion. 

4. While the age determined for M92 ap- 
pears to be too large, within the error bars, 
it is consistent with the recent WMAP results 
of the age of the un iverse of 13.7 ± 0.2 Gyr 



(jSpergel et al. 1120031 ). Accounting for the er- 
ror bars, it is also consistent within 1 — a with 
the mean age of the oldest globular clusters: 
12.6 ± 1.2 Gyr, based on the luminosity of the 
turn-off in Krauss fc Chabover ( 20031 ). 



Based on the observations made with the 
NASA/ESA Hubble Space Telescope, ob- 
tained from the Data Archive at the Space 
Telescope Science Institute, which is oper- 
ated by the Association of Universities for 
Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations 
are associated with program GO-9453. 
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The Luminosity Functions 



Table 1 

The Ground-Based Luminosity Functions 
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Table 1 — Continued 
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Table 2 

The HST Luminosity Functions 
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Table 2 — Continued 
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,24 





.91 


3.48 


20.1 


2419. 


.77 


72 


.77 


0, 


,89 


3 


.38 


3399 


,52 


90 


,58 





.81 


3.53 


20.3 


2502. 


.94 


74 


,27 


0, 


,89 


3 


.40 


3081 


,61 


83. 


,45 





.84 


3.49 


20.5 


2475. 


,90 


73 


.40 


0, 


.90 


3 


.39 


2933 


,70 


79. 


,82 





.84 


3.47 


20.7 


2496. 


,91 


74 


,58 


0, 


,87 


3 


.40 


2395 


,66 


68. 


,53 





.86 


3.38 


20.9 


2730 


,19 


80 


.55 


0. 


,82 


3 


.44 


2121 


,00 


62. 


.69 





.86 


3.33 


21.1 


2572. 


.96 


77 


.24 


0, 


,83 


3 


.41 


2191 


.36 


66. 


,77 





.77 


3.34 
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Table 3 

The Combined Luminosity Functions 



Mag Ny cry compv log(Ny) Ni u\ compi log(Ni) 



10.9 






11.1 






11.3 






11.5 






11.7 






11.9 






12.1 


8.05 


8.05 


12.3 


7.28 


4.21 


12.5 


7.91 


3.96 


12.7 






12.9 


2.49 


2.49 


13.1 


2.17 


1.54 


13.3 


3.22 


1.86 


13.5 


11.05 


4.94 


13.7 


15.83 


5.01 


13.9 


17.79 


5.63 


14.1 


17.43 


4.36 


14.3 


16.27 


4.35 


14.5 


34.18 


6.35 


14.7 


31.71 


5.99 


14.9 


30.05 


5.68 


15.1 


23.55 


4.91 


15.3 


28.51 


5.39 


15.5 


47.83 


6.98 


15.7 


37.47 


6.16 


15.9 


56.06 


7.56 


16.1 


60.31 


7.85 


16.3 


82.59 


9.18 


16.5 


104.43 


10.34 


16.7 


85.20 


9.30 


16.9 


126.53 


11.36 


17.1 


144.65 


12.23 


17.3 


158.51 


12.73 


17.5 


215.00 


14.87 


17.7 


295.86 


17.49 


17.9 


497.36 


22.82 


18.1 


793.09 


28.96 


18.3 


1009.64 


32.71 



1.50 
9.51 
12.23 







2.49 






2.17 


0.12 


0.91 


1.18 


0.41 


0.86 


10.19 


0.51 


0.90 


11.86 






5.58 


0.40 


0.40 


16.84 


0.92 


0.34 


20.97 


0.93 


0.51 


10.18 


0.45 


1.04 


25.50 


0.63 


1.20 


35.45 


0.56 


1.25 


27.41 


0.92 


1.24 


16.49 


0.86 


1.21 


26.62 


0.85 


1.53 


38.66 


0.88 


1.50 


35.58 


0.93 


1.48 


34.38 


0.98 


1.37 


59.16 


0.98 


1.45 


61.34 


0.98 


1.68 


85.56 


0.99 


1.57 


86.16 


0.98 


1.75 


88.11 


0.98 


1.78 


119.41 


0.98 


1.92 


137.33 


0.98 


2.02 


148.79 


0.99 


1.93 


181.88 


0.98 


2.10 


186.55 


0.97 


2.16 


278.08 


0.98 


2.20 


389.58 


0.97 


2.33 


597.74 


0.97 


2.47 


798.94 


0.96 


2.70 


1154.05 


0.95 


2.90 


1506.07 


0.94 


3.00 


1907.65 



1.50 


0.67 


0.18 


5.49 


0.32 


0.98 


6.12 


0.33 


1.09 


2.49 


0.40 


0.40 


1.54 


0.92 


0.34 


1.18 


0.85 


0.07 


4.56 


0.49 


1.01 


4.84 


0.51 


1.07 


2.50 


0.90 


0.75 


5.61 


0.53 


1.23 


4.81 


0.91 


1.32 


3.39 


0.88 


1.01 


5.32 


0.90 


1.41 


6.58 


0.82 


1.55 


5.48 


0.91 


1.44 


4.12 


0.97 


1.22 


5.22 


0.98 


1.43 


6.27 


0.98 


1.59 


6.01 


0.98 


1.55 


5.90 


0.99 


1.54 


7.77 


0.98 


1.77 


7.92 


0.98 


1.79 


9.34 


0.98 


1.93 


9.40 


0.98 


1.94 


9.45 


0.99 


1.95 


11.04 


0.98 


2.08 


11.91 


0.97 


2.14 


12.36 


0.98 


2.17 


13.67 


0.97 


2.26 


13.79 


0.98 


2.27 


17.02 


0.96 


2.44 


20.25 


0.95 


2.59 


25.10 


0.95 


2.78 


29.13 


0.94 


2.90 


34.91 


0.95 


3.06 


40.25 


0.93 


3.18 


45.36 


0.93 


3.28 
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Table 3 — Continued 



Mag Ny cry compv log(Ny) Ni u\ compi log(Ni) 



18.5 


1305. 


.95 


37 


.33 


0. 


.94 


3 


.12 


2190.85 


49.24 


0.90 


3.34 


18.7 


1585. 


.66 


41 


.33 


0. 


.93 


3 


.20 


2428.96 


51.88 


0.90 


3.39 


18.9 


1908. 


.09 


45 


.23 


0. 


93 


3 


.28 


2916.14 


56.89 


0.90 


3.46 


19.1 


2057. 


,72 


47 


.83 


0. 


.90 


3 


.31 


3410.82 


62.72 


0.87 


3.53 


19.3 


2295. 


,62 


50 


.40 


0, 


.90 


3 


.36 


3909.42 


67.44 


0.86 


3.59 


19.5 


2671. 


.86 


54 


.73 


0. 


.89 


3 


.43 


4688.35 


75.42 


0.82 


3.67 


19.7 


3100 


,09 


59 


.62 


0, 


,87 


3 


.49 


4890.62 


76.79 


0.83 


3.69 


19.9 


3331 


,05 


61 


.73 


0. 


.87 


3 


.52 


5108.86 


77.97 


0.84 


3.71 


20.1 


3887. 


.75 


68 


.44 


0. 


,83 


3 


.59 


5669.78 


85.82 


0.77 


3.75 


20.3 


4169. 


.43 


70 


.99 


0. 


,83 


3 


.62 


5718.65 


88.78 


0.73 


3.76 


20.5 


4081. 


.43 


69 


.98 


0. 


.83 


3 


.61 










20.7 


4158. 


.33 


70 


.97 


0. 


,83 


3 


.62 










20.9 


4402. 


.63 


74 


,05 


0. 


.80 


3 


.64 










21.1 


4744. 


,25 


81. 


.59 


0, 


,71 


3 


,68 
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